Dynamics of Bianchi IX universe with massive 

scalar field 



A. V. Toporensky 1 * and V. O. Ustiansky 1 ^ 

1 Sternberg Astronomical Institute, Moscow University, Moscow, 119899, Russia 



Abstract 

The dynamics of the Bianchi IX cosmological model with mini- 
mally coupled massive real scalar field is studied. The possibility of 
non-singular transition from contraction to expansion is shown. A 
set of initial conditions that lead to non-singular solutions is studied 
numerically. 
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1 Introduction 



During last years the dynamics of closed FRW universe with massive mini- 
mally coupled real scalar field has become a subject of investigations. Being 
one of the simplest cosmological models it has a peculiar property. Namely 
the universe of this type can make a non-singular transition from contraction 
to expansion (in other worlds a scale factor can posses a local minimum). 
Such a transition is often referred to as 'bounce'. 

It turns out that a universe in principle can make an arbitrary number of 
bounces. Page was the first who conjectured the presence of a fractal set 
of infinitely bouncing solutions, later Cornish and Shellard |2j proved that 
this system is really chaotic. In later works , fill the chaotic properties 

of this and more complicated models were studied. 

Numerical investigations have shown that the initial conditions space has 
a quasi-periodical structure : there are narrow zones which lead to bouncing 
solutions separated by wide zones which lead to singular solutions. These 
investigations also have shown that in former zones the solutions exhibits a 
sensible dependence of initial conditions - a property peculiar to the chaotic 
dynamical systems. Moreover, starting from one of such zones a solution may 
pass after bounce through another zone of this type (or through the same 
one) and suffer a new bounce. After each bounce the situation is similar 
to the previous one. This lead to the existence of a fractal set of unstable 
infinitely bouncing solutions. In the dynamical chaos theory such set is called 
as a strange repellor - a structure which can be found with in many chaotic 
dynamical systems without dissipation. The presence of this fractal structure 



of solutions corresponds to the fractal structure of sets of initial conditions 
leading to bouncing solutions @. 

But the Universe at the earliest stages of its evolution have not necessary 
had a metric close to FRW one. So for studying the Universe at this time 
one often invoke more general metrics. The simplest non-FRW metric with 
positive curvature is the Kantowski-Sachs one. However this case does not 
include the closed FRW metric, and moreover there are no bouncing solutions 
in this case (see appendix). More complicated metric is Bianchi IX model. 
This is the only model that includes closed FRW metric ||. So in this paper 
we investigate the Bianchi IX universe with massive scalar field. 

The Bianchi IX universe is well known due to its mixmaster approach 
to singularity ||. It worth to notice that mixmaster chaos 0, |J principally 
differs from the studying one as in the case of mixmaster it is the evolution 
of anisotropic variables which exhibit chaotic properties while the overall 
volume of the universe decrease monotonically. In the considered case it 
is the evolution of the total volume of the universe which exhibit chaotic 
properties. 

2 Equations of motion 

Let us consider a Bianchi IX universe with minimally coupled massive scalar 
field with the action as follows (we use units c— h — 1, G — I/871"): 

S = J (-f + \w* - d*x . (1) 
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The metrics has the following form : 

ds 2 = dt 2 - r) a p{t)uj a ujP , 

where 



(2) 



a(t) 2 \ 
Vcfi(t) = | b{t) 2 , (3) 
c(t) 2 J 

and u a are differential 1-forms invariant under SO (3) transformations. 

Assuming ip to be a function of time variable only one can write the 
Einstein's field equations in following form : 



(abc) 1 , 4 , 2 o\2i m2 ,„2 

( Q M i j; r^4 / 2 2\2l m2 , 2 

a&c 2a 2 b 2 c 2 I 

a b c , 

<p+ [ - + 7 + - ) tp + m z ip = 
a o c 



(4) 

(5) 
(6) 

(7) 



with the first integral of motion : 
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where 

In our investigation it is more convenient to work with ADM variables 
{Q(r), /3 + (r), /3_(r)} which are functions of a new time variable r : dt — 
m 5 e 3n dr . Relations between ADM variables and the old ones are : 
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It is easily seen from ([]) that variable Q describes the volume of the universe 
while f3± describes it's anisotropy. 

The considered dynamical system admits Hamiltonian formalism. In 
terms of ADM variables the Hamiltonian has the following form : 
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h = -(-pI + pI+pI + -pi) + -e 4n y(/3 + ,/3_) + -<A 6n , (io) 



where 



V(/3+,P-.) = \{e~ m +2e 4/3 + (cosh Ay/30- - 1) -4e~ 2/3 + cosh2v / 3/3_) . (11) 

The Hamiltonian equations leads to the following : 

n" = i e 4 V(/3 + ,/3_) + ^V n , (12) 

f3'l = \ e m \e~ 8 ^ - e 4/3+ (cosh4v / 3 / 9- - 1) - e' 2 ^ cosh 2^.1 , (13) 

^/ = L-2^ +siim2 V3/3_ -e 4/3 +sinh V3/3-1 , (14) 

3 - 

= -<^ e 6C (15) 



and the constraint equation H = takes the form : 

- Q' 2 + /T 2 + + ly?' 2 + ^W(/3 + , /3_) + = , (16) 

where ' = . 

To make further investigations we need to generalize the notion of bounce 
to Bianchi IX case. Being defined as a non-singular transition from contrac- 
tion to expansion of the universe bounces thus corresponds to local minima 
of Q. That is why ADM variables is more preferable choice for our purposes. 
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Substituting Q' = into constraint equation ( |TED and combining it with 
the condition Q" > it is easily seen that bounces can occur only if : 

V(P + ,f3.) + ^ 2 e 2n <0, (17) 
V(P + ,P_) + ^ 2 e 2n >0. (18) 

The above inequalities at least implies V((3+,(3J) < which can be satisfied 
only in a narrow region around the isotropic case (j3 + = (3- = 0) where 
V(/3+, suffer it's minimum. This region is represented on fig. 1 where 
the levels V(/3+,/3-) = const are shown for the values of const = — .99(inner 
curve), const = —.90 , const = —.75 , const = —.50 and const = 0(bold 
lines). The potential has a minimum in point /3 + = (3- = (isotropic case) 
and V(0, 0) = —1. It also follows from ( |T6| ) that if a bounce occur the possible 
values of f3'± must satisfy criterion : (3+ + (3® < 1 as V((3+,(3-) > —1 for all 
values of f3± . 

Before starting the next section it should be mentioned that Belinsky and 
Khalatnikov [|TJ were the first who studied the Bianchi IX model whith scalar 
field. Thay considered only massless case and found that the presence of the 
scalar field of this type destroys the Mixmaster oscillatory regime when the 
universe collapses to singularity. This results holds in massive case as an 
asymptotic regime. 

3 Numerical results 

The equations of motion QI^D-flOp can be solved numerically. In order to 
obtain numerical solution one need to fix eight initial conditions. The con- 
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Figure 1: A levels of V(/3+,/3_) = const. 



straint equation allows to express one of them in terms of the others (in our 
work we determined 0(0) in this way). Moreover as the considered system 
has at least one maximal expansion point (i.e a point where Q suffer (local) 
maximum), one can fix fi'(O) = without loss of generality. 

For simplicity in our calculations we always fixed ip(0) = . So we 
studied the dependence of numerical solution from five initial conditions : 
fi(0), /?+(()), /3_(0), P' + (0), P'_(0) . 

The equations of motion have been solved numerically using Bulirsch- 

)-fll5|) requires evaluations of transcen- 



Stoer method |10| . As equations (|T 
dental functions while equations (^)- (^) requires only arithmetical opera- 
tions we performed calculations in variables {a, b, c} in order to speed-up 
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them while the initial conditions and results were expressed in terms of 
{Q, /3 + , j3-} . The constraint equation @ has been used to check the ac- 
curacy. 

First of all we tested our routines on already known fridmanian case. 
On fig. 2 we presented a first four closest to singularity simplest symmetric 
periodical solutions in a projection onto (<£>, Q) plane. Term 'simplest' means 
that all of the maximal expansion points for given solution coincide in this 
plane and symmetric means that these points lies on axis ip = 0. The values 
of Q in the above mentioned points are : Q± = 0.160832 (short dashed 
line), Q 2 = 0.760135 (long dashed line), Q 3 = 1.141779 (dot-dashed line), 

= 1.418301 (solid line). Near each of these values fl{ there are narrow 
intervals which lead to bouncing solutions. 
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Figure 2: Examples of periodical solutions in Fridmanian case. 

As Q'(0) = the initial values /3±(0) and (3±(0) have similar restrictions as 
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in the case of bounce with the only exception that (|TSD is no longer required. 
Numerical investigations have shown the presence of bouncing solutions only 
in the case of sufficiently weak anisotropy near the corresponding fridma- 
nian bouncing solutions. These investigations have also shown that for some 
ranges of initial values of Q, f3± and (3'± solutions exhibits sensible dependence 
of initial conditions. 

To illustrate the structure of the sets of initial conditions which leads to 
bouncing solutions we calculated a 2d slices through the initial conditions 
space. For given value of Q(0) each slice represents a 500 x 500 grid in a 
plane of two of the anisotropy variables while the rest of them equals to zero. 
For each knot of the grid it has been calculated a numerical solution on a 
sufficiently large interval of t G [0, 50] . If a maximal value of Q on this 
solution after bounce fl ma x is grater then f2 4 we encoded corresponding dot 
in a black color. If f2 4 > Q max > f2 3 we encoded it in a dim gray color. 
If ^3 > Qmax > ^2 we encoded it in a gray color and if ^2 > &>max > ^i 
we used a light gray color. We considered only solutions which satisfies 
criterion VL max > Q\ as they are of the most interest. We do not consider 
the trajectories with the next maximum expantion point lying closer to the 
singularity than the first bounce interval, becouse in this case the trajectory 
have only a zigzag, which can not change the fate of the trajectory falling 
into the singularity. 

We calculated such slices for subsequent values of Q with sufficiently small 
steps so that they can be viewed as a 'movie'. Such movies can be found 
at "http: / / www.xray.sai. msu.su/~ustiansk/sciwork/gr-qc/001. html" . Exam- 
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pies of these slices are shown in fig. 3-4. 

Numerical calculations have shown that there are narrow intervals of ini- 
tial values Q that lead to bouncing solutions separated by wide intervals that 
lead to singular solutions in full analogy with the fridmanian case. For the 
closest to singularity interval (an interval which contain Qi) in a 
plane one can observe the picture that follows. For the values f2 < 0.1549 
there are no bouncing solutions. With the increase of Q bouncing solutions 
appears in a small region around isotropic case which grow with the growth 
of up to the value f2 ps 0.1567 when a fraction of considered anisotropic 
bouncing solutions begins to decrease until they vanishes at Q m 0.1608 . 

In the case of (f3 + , f3' + ) and (f3~, f3'_) slices the general picture slightly 
differs. With Q increasing from the value f2 ps 0.1525 anisotropic bouncing 
solutions appears and invoke isotropic case at already mentioned value Q ~ 
0.1550. Thus in considered case the interval of Q which allows bouncing 
solutions becomes slightly wider than in fridmanian one. 

For the interval which contain f2 2 the picture is following. For both 
studied planes and (f3-,f3'_) bouncing solutions appears when VL ps 

0.7593 in a small region around isotropic case. With the growth of Q this 
region grows until Q ~ 0.7603 . After that the value of Vt max becomes smaller 
than Qi for initially anysotropy close to zero so such solutions were discarded 
and thus the region have a 'ring' shape. For the values VL ps 0.762 this regions 
vanishes. So in considered case the interval of Q is wider than in fridmanian 
one. The similar picture is valid for the interval that contain Q 3 . 

On several slices one can easily see a narrow 'rings' (one of them is pointed 
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by two arrows in fig. 4) inside the region of initial conditions that lead to 
bounces. The existence of these 'rings' is caused by the fact that solutions 
which correspond to these 'rings' suffers bounce and the next point of max- 
imal expansion also occurs inside one of the regions that leads to bounce. 
So such solutions suffer at least two bounces. These 'rings' forms a fractal 
structure. The example of such structure is shown in fig. 5. 

For the splits in (f3 + , (3 J) it is possible to determine the mesure of bouncing 
solutions in a simple way if we start from </?(0) = 0. In this case the measure 
of all possible initial values of (3± , being defined as the area bounded by the 
bold lines in fig.l , equals ~ 0.48. The maximal fractions of initial values 
of f3± that lead to acceptable bouncing solutions to all posible initial values 
for first three subsequent intervales of f2 are approximately 0.0026, 0.00041, 
0.00012. 

4 Conclusions 

In this paper we studed the Bianki IX cosmological model with the massive 
scalar field. In particular the influence of initial shear on the possibility of 
transition from contraction to expansion was investigated. We show that only 
a very restricted set of initial conditions with shear can led to the bouncing 
trajectories. Their mesure with respect to all physically admissible initial 
conditions for the first 3 zones are restricted by 2.6 x 10~ 3 , 4.1 x 10~ 4 and 
1.2 x 10~ 4 respectively. Corresponding areas of the initial conditions were 
found numerically also for other cross-sections of the initial conditions space. 
In all cases, bounces are possible for a ruther narrow zones near corresponding 
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friedmannian trajectory. 

On the other hand, the structure of chaos is similar to known in the 
Friedmann case: there are narrow intervals of initial values Q that lead to 
bouncing solutions separated by wide intervals that lead to singular solu- 
tions. The width of the bounce intervals may be enlarged up to 3 times in 
comparition with the Friedmann case for the several closest to the singularity 
intervals, but the qualitative picture is the same one. 

Appendix. No chaos in Kantowski— Sachs uni- 
verse with a massive scalar field 

The equations of motion for the Kantowski-Sachs universe filled by a scalar 
field are 

b ab <b 2 . . 

b ab a 6 2 Tr/ ,. . . 

l + H + a + T = V ^ < 2 °» 

4>+(- + 2\)4> + V'{4>) = 0, (21) 
a 



with the constraint 



hb b 2 1 4> 2 t r / i \ 

2 s + ^ + ^ = y + W) (22) 



where a, b - the scale factors, V((f>) - a potential of the scalar field, and the 
prime indicates the derivative with respect to 0. 

It is easy to see from the Eq.(19) that in all the possible turning points 
with respect to b (i.e. the points in which 6 = 0), the corresponding acceler- 
ation 

b = -b<j) 2 /2 (23) 
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is always nonpositive, so we can have only the points of maximal expansion 
for the scale factor b. 

Combining Eqs.(19 - 20) we receive, on the other side, that the acceler- 
ation in the turning points with respect to a is 

a = aV(<f>), (24) 

which is nonnegative for an arbitrary nonnegative scalar field potential V(<p), 
and we have only the points of minimal contraction for the scale factor a. 

As a result, we can claim from Eqs.(23 - 24) that the type of chaos asso- 
ciated with the oscillation of the scale factors is absent in Kantowski-Sachs 
universe filled by the scalar field with an arbitrary nonnegative potential. 
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Figure captions 

Fig. 3: A slice through initial condition space with fixed f2 = 0.1573 which 
represents a 500 x 500 grid in (j3 + , f3S) plane. For each knot of the grid it has 
been calculated a numerical solution. If a maximal value of Q on this solution 
Qmax > ^4 we encoded corresponding dot in a black color. If f2 4 > Vt max > Q 3 
we encoded it in a dim gray color. If Q 3 > Vt max > Q 2 we encoded it in a 
gray color and if Q 2 > ^max > &i we used a light gray color. 

Fig. 4: A slice analogous to previous one for the value f2 = 0.7604. One of 
thing rings which represents second order structure is pointed by two arrows. 
These rings are formed by trajectories which suffer two bounces. 

Fig. 5: Example of fractal structure in plane with the value of 

Q = 0.7604. This is a magnified region pointed by two arrows on previous 
figure. Wide regions represents a second order structure. One can easily see 
a third order structure as thing lines. 
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